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Inverse-Scattering Theory and the Density Perturbations from Inflation 
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Inflationary cosmology provides a successful paradigm for solving several problems, notably the 
generation of density perturbations which seed the formation of observed cosmic structure. We show 
how to use inverse scattering theory as the basis for the inflationary reconstruction program, the goal 
of which is to gain information about the physics which drives inflation. Inverse scattering theory 
provides an effective and well-motivated procedure, having a sound mathematical basis and being 
of sufficient generality that it can be considered the foundation for a non-parametric reconstruction 
program. We show how simple properties of the power spectrum translate directly into statements 
about the evolution of the background geometry during inflation. 

PACS numbers: 98.80.Cq 



A wealth of data drives much of the current activity in 
cosmology. Amongst these data, perhaps the cleanest sig- 
nal from the primordial Universe is the power spectrum 
of density fluctuations, which is obtained from observa- 
tions of fluctuations in the microwave background sky 
and from observations of the large-scale clustering of mat- 
ter 0]. Especially in the case of microwave background 
observations, the observed fluctuations give direct infor- 
mation on the "initial conditions" for the density pertur- 
bations which manifest in the gravitational clustering of 
matter in the Universe. 

In an inflationary universe, these initial density per- 
turbations arise as a relic of quantum fluctuations from 
the very earliest times. The inflationary evolution causes 
these modes to grow in amplitude, with a time-evolution 
that is determined completely by the evolution of the 
space-time geometry, encoded in the scale- factor a{t). 
This evolution imprints a signature of the inflationary 
evolution on the spectrum of modes. Therefore, struc- 
ture in the mode spectrum can be directly related to the 
evolution of a{t) and hence to the nature of the stress- 
energy which drives inflation. Because the nature of in- 
flation is quite mysterious, information of this sort is very 
important for cosmology. 

Consider now the standard picture for quantum fluc- 
tuations and the generation of structure during inflation. 
The equation describing the evolution of gauge-invariant 
metric fluctuations can be taken in the form [3l| 



< + [e - qir,)] <t>k = 0, 



(1) 



where primes indicate derivatives with respect to con- 
formal time 7], and (7(77) is an effective potential. The 
spatial dependence of the functions has been expanded 
in modes specified by three-momentum fc, so that 4>k is a 
function of 77 only. For scalar modes we have 9(77) = z" / z, 
where zijf) is a function determined by the evolving back- 
ground; under quite general circumstances it is given by 



z{ri) 



airj) 



where h{ri) = a'{rf)/a(ri) and Cs{ri) is the "adiabatic 
sound speed" . The curvature term K. rapidly becomes 
negligible as inflation proceeds and will be ignored for 
the remainder of this Letter. When inflation is driven 
by a scalar field, = 1. For tensor modes the effective 
potential is given by qriv) — o." {v) / ^{"H) ■ the infinite 
past, the mode population is that of a precisely specified 
vacuum state, which gives the condition 



cj)k ~ {2k)- 



The argument of the exponential is determined by the 
positive-energy condition in the deep past, and the nor- 
malization is determined by the usual free-field normal- 
ization. The constant phase factor 7 is a matter of con- 
vention and is fixed later in order to simplify resulting 
expressions. 

A given fc-mode evolves to a good approximation as a 
free mode until it begins to feel the effective potential, for 
rj such that ~ liv)- When it crosses into the forbidden 
region (crosses outside the Hubble length), it begins to 
grow in amplitude. Deep within the forbidden region, it 
can be easily shown that the scalar modes behave like 
(f'k ^ ^kzii]), where Af^ is a function of k only. The 
structure of q{r]) is imprinted on this function, which is 
directly related to the power spectrum of the fluctuations 



P{k)^^\Ak\ 



,{-q)h{7]) 



[h'{r^)-h'{r^) + K] 
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This power spectrum is in turn directly related to the 
spectrum of density fluctuations or curvature fluctua- 
tions. 

The goal of inflationary reconstruction is to determine 
as much information as possible about the stress-energy 
which drove the inflationary expansion. The reconstruc- 
tion program splits naturally into two clearly defined 
tasks. The first task is to solve the linear problem, in- 
verting from the observed power spectrum to the func- 
tion q{rj). The second task is to use the obtained q{rj) 
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to constrain the geometry and therefore to constrain the 
stress-energy tensor which drives inflation. In this Letter 
we show how to use inverse scattering theory to complete 
the first task; we then show how this information can be 
used in the second task, to constrain the physics of the 
inflationary epoch. 

The basic problem of calculating density perturbations 
is most naturally posed as a problem in scattering theory 
for the wave equation Q ■ In the following we can restrict 
ourselves to the case of an effective radial Schrodinger 
equation for scattering from a central potential. In par- 
tial wave I this equation is 



0"(fc,r) 



1) 



V{t) 



(j){k,r) = 0. 



In order to avoid confusion in applying the formalism to 
the wave equations given above, we identify r — —rj and 
work solely in terms of the variable r. 

The reason to concentrate on the partial wave equation 
at fixed £ is the following. It turns out to be most natural 
to split the potential in the form 



q{r) 



where 9 is a judiciously chosen constant. We identify 
£(^-1-1) = d'^ — 1/A and the remaining term as the central 
potential V{r). In this way we collect the most singu- 
lar part of the potential into the effective "centrifugal" 
term. The remainder of the discussion explains precisely 
how this process is carried out and what it means. Note 
also that there are no bound states in the problem, inde- 
pendent of the sign of V{r), because of the positivity of 

The wave solution is purely an incoming wave for 
r — > oo, which follows from the positive-energy vac- 
uum condition in the deep past. Solutions which obey 
this condition are important in scattering theory and are 
called Jost solutions. Precisely, a Jost solution is the 
unique solution of the wave equation satisfying the con- 
dition 4>{k,r) ^ exp{ikr + i7r£/2), r oo. Note that a 
Jost solution must become singular deep inside the for- 
bidden region, as a consequence of the lack of reflection in 
its definition. In our case the Jost-like solution is actually 
the physically interesting solution, describing a positive- 
energy quantum mode from the deep past. 

The behaviour of a Jost solution deep inside the for- 
bidden region is described by a function of /c. 



FAk) ^ lim 



^r(i) 



r(^+i) V 2 ' 



where fg{k,r) is the Jost solution for k. The function 
F£(fc) is the so-called Jost function, and it encodes es- 
sentially all the information about the scattering theory. 
Given our setup of the problem, we have z ~ co{ro/rY for 
r — >• 0, where cq is a fc-independent constant. Therefore 



a direct relation holds between the Jost function Fe{k) 
and the basic function of interest Ak, 



Ak = 



Cora) 



f+1/2 



Fe{k). 



At this point it is possible to invoke the full machinery 
of inverse scattering theory, which allows us to recon- 
struct the potential V{r) from the Jost function, by solv- 
ing a linear integral equation 0| . Rather than solve this 
numerical problem, in this Letter we instead show how 
the Jost function encodes information about the asymp- 
totic behaviour of V(r), and therefore the function z{r), 
in a manner which is directly applicable to the cosmolog- 
ical inversion problem. 

To be precise, suppose that we write z(r) ~ 
(?'o/f')^c(r), where c(r) is assumed to be regular and to 
satisfy c{oo) = 1. The potential is then given by 



Vir) 



z 
z 



^(^ + 1) 



2^c^ 
r c 



By direct calculation, we can show that Fi{0) — c(0). 
This is a consequence of the special form for V{r). A 
more difficult result expresses the behaviour of Fe{k) for 
fc — > in terms of the large r behaviour of V{r). To this 
end we define the length-scale ri and the exponent S by 



Vir) 



( r ) 



The small k behaviour for £ = has been calculated by 
Klaus 0. Extending Klaus' calculation to general i, we 
obtain the fc — > asymptotic 



Fejk) 
Fi{Q) 



= 1 



T(5/2 



r(^-i- 1)2^- 



(2) 



where A{£, S) = 6) + A2{£, 5) with 



= 2 



'-1-5' 



Ai{£,5) 
A2{£,5) = V(£-5) 



i)r(i 



T{\)T{2 + £ + 5) 



This particular calculation assumes < (5 < 1, which 
turns out to be the generic case for inflationary models; 
other cases can be handled in a similar manner. Given 
this result, we see that knowledge of the Jost function for 
k near zero translates into knowledge of V{r) and z{r) for 
large r. This result is all we will need for the remainder 
of the discussion. 

Before proceeding let us examine the meaning of the 
various ingredients introduced thus far. First consider 
the form chosen for the potential. In power-law inflation, 
where the scale- factor is assumed to grow like a{t) ~ , 
the function v{r) is a constant, v = 3/2-1- l/{p~ 1). The 
fact that the power-law case corresponds precisely to a 
free wave equation is essentially a statement of the exact 



scale-invariance of the wave equation in that case. The 
usual spectral index and the partial wave index are re- 
lated by n = 3 — 2£. The Harrison- Zeldovich case with 
n — 1 corresponds to partial wave i = 1. When £ is dif- 
ferent from zero, the power-spectrum has a pure power- 
law scaling form with an anomalous scaling exponent. It 
is interesting to note that the Harrison- Zeldovich case 
corresponds to an anomalous exponent, although from 
another point of view it is "natural" . 

Only deviations from scale-invariance in the wave 
equation can give rise to deviations from power-law be- 
haviour in the spectrum of density perturbations and to 
a nontrivial Jost function. Clearly the Jost function is 
identically one in the free case. Note that the parame- 
ter £ need not be an integer; nothing in the analysis is 
affected by this. 

Next consider the meaning of the limit ?' ^ which 
defines the Jost function and the spectrum of density 
perturbations. Obviously the detailed structure of the 
potential in a neighbourhood of the origin depends on 
the way in which inflation ends; if the details of this pro- 
cess become important, then we lose the ability to make 
robust predictions about the spectrum of density pertur- 
bations. However, this is simply a statement defining 
the range of k modes of interest. In practice we observe 
a range of modes, which we assume were "well-inflated", 
so that they are not sensitive to such details. A typical 
complementarity relation holds, so that this range of k 
modes corresponds to a range of r. In this spirit, the 
range of k modes of interest corresponds to small k and 
therefore to large r; these are the modes which spent a 
significant amount of time outside the horizon. 

Combining the expressions given above, we obtain a 
representation for the power spectrum. 



10" 



P{k) 



\r{£- 



1)|2 



47r3 



-21 



Fi{k) 



Fi{0) 



(3) 



This form for the power spectrum illustrates several ideas 
in the scattering theory approach. The partial wave in- 
dex £ is directly related to the anomalous scaling expo- 
nent in the leading power-law behaviour, and the momen- 
tum scale for the anomalous part of the scaling relation 
is set by r^^ . These observations completely specify the 
meaning of the parameters £ and rg in this formalism. 
All deviations from scaling are contained in the shape 
function, \Fi{k)/Fe{0)\. 

The usual form of the power-spectrum given in the 
literature is non-dimensional, being a power-spectrum for 
a non-dimcnsionalized curvature scalar TZ jy]. Relating 
the definition of TZ to the objects given here, we find 

Pizik) = 2TiGP{k). 

This yields an object with no engineering dimensions. 
Note that we have set Ti and c to unity everywhere, so 
that G has units of length-squared. At this point it is 
possible to interpret a large class of power-spectra in a 
model-independent way. Next we illustrate this approach 
by example. 
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FIG. 1: Power spectrum computed for a <j) inflation model. 
The dashed line indicates the small k behaviour, a pure power 
law. The fitting procedure is discussed in the text. 



As an example, consider the power spectrum in Fig. 
n This power spectrum was computed numerically from 
an inflation model with a (j)^ inflaton potential. Com- 
putations of power spectra for other single-field mod- 
els show that the properties of this power-spectrum are 
generic for a large class of such models and it serves as a 
good illustration. The behaviour as k approaches zero 
allows us to directly determine the parameters £ and 
tq. The global behaviour is fit to the fc — > asymp- 
totic form given by the combination of Eqns. Q and 
(PJ. For fitting, one must specify a range of k val- 
ues in which the asymptotic result is assumed to hold. 
Using the full range of data for the power spectrum, 
which is displayed in Fig. we obtain £ ~ 1.02496, 
5 ~ 0.103, and ri ~ 6.43 x IQ-^^ /j-i]y[p(,sec/km. Us- 
ing a significantly more restricted set of the data, with 
k < 0.0001 /iMpc~\ we obtain £ ~ 1.02603, S ~ 0.118, 
and ri ~ 1.39 x 10~^^ /i~^Mpc sec/km. The large diflFer- 
ence in the determinations for ri reflects the fact that 
the exponent 5 is small; the actual coefficients of the cor- 
responding potentials differ only by about 10%. 

With this information we reconstruct the large r be- 
haviour for V{r) and therefore z(r). In Fig. |2]we com- 
pare the reconstructed function z" {r) / z{r) — r'^V{r)-\- 
t{£ + \) to the input which gave the power spectrum of 
Fig. n Note that both reconstructions agree well with 
the input and with each other for large r, as they should. 
A small additive error results from lack of precision in 
the determination of £, which determines the pedestal for 
each of the functions. The input function rolls over for 
small r, which indicates the end of infiation; of course, 
this rollover is not captured by the asymptotic recon- 
struction. Note further that the amplitude of the power 
spectrum has no effect on these results; the amplitude 
simply determines the scale tq. In this case tq is roughly 
5 X 10^ G^/^. The fact that tq is much larger than the 
Planck length reflects the fine-tuning necessary to obtain 
a viable inflationary model. 
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FIG. 2: Comparison of reconstructed scattering potential to 
the input. The input was tabulated during computation of 
the power spectrum. 



To close the argument, we relate the large r asymptotic 
for V{r) to the early time behaviour of the Hubble ex- 
pansion. If we assume a past asymptotic for the Hubble 
expansion in the form 



H{t) =pt-^ 



t 0, 



then a short calculation shows that 



Vir) 



-(-) 

j,2 \ J. J 



where ri can be written in terms of t^, though the relation 
is not needed here. Applying this result to the example 
at hand, we obtain h ~ 4.2 — 4.5. The large power is 
consistent with the idea that the evolution of the scale 
factor for our (fP' model is very close to power-law inflation 
for most of the inflationary epoch; this is why the power 
spectrum is close to a simple scaling law. 

Thus far we have discussed the power-spectrum for 
scalar fluctuations. We now briefly show how to relate 
the scalar and tensor spectral indices using the behaviour 
of the equation of state. 

The stress-energy driving inflation can be parameter- 
ized by an equation of state function w(r]), defined by 
p = wp. Recall that one of the important predictions 
of inflation is a relation between the spectral indices for 
scalar and tensor fluctuation modes. During inflation the 
equation of state is close to w = — 1. As inflation grad- 
ually ends, the equation of state drifts away from —1. 



Assume a behaviour of the form 



,11 



We consider the simplest case for 2(77), equivalent to as- 
suming that inflation is driven by a scalar-field. In this 
case it is easy to show that the following relation holds: 

a' z' I w' 



a z IX-V w 
Using this relation, a short calculation allows us to iden- 
tify the partial wave index for tensor perturbations as 
It — I ~ where I is the value for scalar pertur- 

bations. Therefore we obtain a statement of the scalar- 
tensor relation in terms of the parameter m, 



riT 



ns - I 



These spectral indices are defined as outlined in the 
above, in terms of the small k behaviour of the power 
spectra, where the Jost function contributions are negli- 
gible. For power-law inflation, i = £t, and we obtain the 
exact result nx = ng — 1. 

The procedure we have described is essentially non- 
parametric. No particular model for the stress-energy 
need be assumed, and geometric information is extracted 
directly from the power spectrum. This analysis does not 
depend on slow-roll parametrizations or assumptions, as 
distinct from previous treatments 0, 01 ■ Equation (0) 
provides a very clear separation of the contributions to 
the power spectrum, and this separation enables the pro- 
cedure. Furthermore, the relations of scattering theory 
provide a very convenient general framework for posing 
the density-perturbation problem in inflationary cosmol- 
ogy 

We have seen that inflationary reconstruction is es- 
sentially a problem in inverse scattering theory. Direct 
relations were used to invert from the power spectrum 
to asymptotic quantities of interest, which describe the 
background evolution through most of the inflationary 
epoch. Although the full machinery of the Gelfand- 
Levitan equation Q was not introduced here, we are 
currently investigating its use for future analyses. For 
example, the Gelfand-Levitan equation can be applied to 
the analysis of any features which do not represent small 
deviations from scaling, such as might occur in multi-field 
scenarios or other dynamically complicated models. 
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